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Abstract 



In a framework of intranuclear cascade (INC) type calculation, we study a mo- 
mentum spectrum in reactions (K~,K + ) at a beam momentum of 1.65 GeV/c. 
INC model calculations are compared with the relativistic impulse approxima- 
tion (RIA) calculations to perform the detailed study of the reaction mechanism. 
We find that the INC model can reproduce the experimental data on various 
targets. Especially, in the low-momentum region, the forward-angle cross sec- 
\Q ! tions of the (K~ ,K + ) reaction on from light to heavy targets are consistently 

explained with the two-step strangeness exchange and production processes with 
various intermediate mesons, and (ft, ao and /o productions and their decay into 
K + K~ . In the two-step processes, inclusion of meson and hyperon resonances is 
O ' found to be essential. 
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1 Introduction 



The study of nuclear systems with strangeness S = -2 is one of the most important 
and hottest subjects in nuclear physics, since these systems give us unique information 
on YY interaction and they might be a doorway to study multi-strangeness systems 
such as the strange matter. Until now, several ways to produce these S — -2 nuclear 
systems have been proposed pj. 

The first one is the S~ absorption reaction at rest on nuclei. In these thirty 
years, there are three reports on the discovery of double-hypernuclei through the S~ 
absorption. These events are found in the emulsion by searching for the sequential weak 
decay processes (double stars) @]. Although the absorption reaction is the most 
effective and the most direct way for the double-hypernucleus hunting, it is necessary 
to produce slow S~ particles, to stop them in the emulsion, and to search for the 
double stars in the emulsion. Thus the information is still too scarce to determine 
YY interaction at present, and further study along this line is necessary. The second 
way to produce S = — 2 nuclear systems is the production of double-hypernuclei at 
the (K~,K + ) reaction vertices on nuclei. If S-nucleus bound states exist and their 
widths are small enough, we can see discrete peaks in K + momentum spectra which 
directly reflect the "EN interaction. In addition to this clear signal, it is expected that 
double-hyperfragments produce after double-hyperon compound nuclear formation and 
evaporation process in the lower K + momentum region [|J]. Actually, KEK-E176 |J 
group have found S — -2 nuclei at the (K~ , K + ) reaction point in the emulsion target, 
although its mass cannot be specified. The (K~,K + ) reaction also provides us with 
H~ particles which can be used as a source of 5~ absorption at rest. Therefore, the 
study of the (K~,K + ) reaction on nuclei is important and urgent to explore S = -2 
nuclear systems. The relativistic heavy-ion collision is considered to be the third way 
to produce S = -2 nuclear systems. In this reaction, more strange nuclei [S < —3) 
or strangelets are expected to be produced. However, at the current stage, there is no 
positive report on the discovery of these nuclei, probably because the production cross 
section is too small to observe. 

Recently, it has become experimentally possible to study double strangeness ex- 
change reactions (K~,K + ) on nuclear targets. For example, small angle (K~,K + ) 
cross sections at px- = 1-65 GeV/c were measured on several targets at KEK-PS by 
Iijima et al. The measured K + momentum spectrum shows a striking structure 
which gives rise to a target mass number dependence (A- dependence): In addition to 
quasifree peaks which correspond to elementary process of K~p — > K + E~, there ap- 
pear a large bump spreading from p K + = 0.35 GeV/c to 1.0 GeV/c, and the yield of 
this bump in the lower momentum region grows much faster (~ A os ) than that of the 
quasifree peak (~ A ' 3 ) as the target mass increases. The comparison of the measured 
spectrum with the DWIA spectrum indicates that the high momentum peak p^+ ~ 
1.1 GeV/c has been understood by quasifree process K~p — > K + E. However, the huge 
lower momentum bump cannot be explained by the quasifree processes only with well- 
known elementary reactions such as K~p — > K + S* ||. Thus, the appearance of this 
lower peak bump has casted a doubt on the belief that quasifree processes dominate 
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in meson induced reactions at around 1 GeV/c region, and at the same time, it gives 
us a hope that we can make more abundant double-hyp ernuclei from the unknown 
mechanism in this large bump. Therefore, it is important to obtain the knowledge for 
the reaction mechanism of the (K~, K + ) reaction on nuclei. 

Recently, one of possible mechanisms is proposed by Gobbi et al. . They have 
claimed that a large part of the lower momentum peak is exhausted by the contribution 
from the heavy- meson production followed by its decay: K~p — > MA, M — > K + K~ , 
where M is the / (975), a (980) or 0(1020) meson. It should be noted that the first 
reaction is a so-called subthreshold particle production. The treatment of the nucleon 
Fermi motion is considered to be important. Since the incident K~ particle has a large 
cross section with a nucleon in nuclear targets, the first-step reaction occurs mainly in 
the surface region where the Fermi momentum is smaller than that in the inner region. 
Although this reduction of the effective Fermi momentum should largely suppress the 
subthreshold particle production, it is not taken into account by Gobbi et al. In fact, 
we can see much less kaon pair creations M — > K + K~ experimentally in the scintilator- 
fiber detector than those expected in their calculations. Therefore, more detailed 
analyses of these reactions are desired before we conclude that the lower momentum 
bump comes from heavy-meson decays. 

Another plausible process is the two-step strangeness exchange and production 
reaction, such as K~p — > M'Y followed by M'N —>■ K + Y, where M' and Y denote 
intermediate mesons (tt, p,rj, . . .) and hyperons (A, S, A*, £*), respectively. In previous 
works |5], |J, the contributions of these two-step processes are estimated to be too small 
to explain the lower-momentum huge bump. However, these estimations limit the 
intermediate mesons and hyperons to pions and ground state hyperons, respectively. 
In addition, they considered only the 12 C target case. As pointed out by Iijima et al. ||, 
other possible reactions may contribute to the K + spectrum significantly. Furthermore, 
the target mass number dependence in the experimental spectrum is well-fitted with 
^o.56±o.02 j n j QW mom entum region, while ^°- 38±0 - 03 i n the high momentum region. 
This scaling property strongly suggests the contribution of multi-step processes ||. 
Therefore, it is necessary to treat multi-step processes correctly to evaluate the K + 
momentum spectrum. 

In this paper, we investigate multi-step effects in the (K~ , K + ) reactions on the 
various targets, i.e., 12 C, 27 A1, 63 Cu, 107 Ag and 208 Pb, using the intranuclear cascade 



(INC) model [|TTJ, [TTJ, [L2[ [y|. One of our aims is to see the target mass number 
dependence on the K + momentum distribution in the (K~,K + ) reaction. Thus as a 
first step, we neglect some other effects such as the influence of the potential effect 
between hadrons. We will also compare the results of the INC spectrum with the one- 
step process spectrum in the relativistic impulse approximation (RIA) model |4lJ to 
study the reaction (K~~,K + ) in detail. 

In order to calculate the contributions of various two-step processes, there are 
two problems to be settled. The first problem is that they involve several elementary 
processes such as pN — > KY, whose cross sections cannot be measured experimentally. 
In this paper, we adopt the generalized Breit-Wigner formula ||, [!|, in which s channel 
dominance is assumed and the interference effects between different baryon resonances 
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are ignored. It has been shown that this formula works well in the energy range 
under consideration, and it reduces theoretical ambiguities. The second problem lies 
in the treatment of complex multi-step processes. For the discussion of multi-step 
processes, cascade type models are considered to be the most reliable ones at the 
current stage: Each elementary two-body collision process is treated explicitly, then it 
is easy to include various multi-step processes. The INC model does quite well for many 



situations. After the first application to high-energy nuclear collisions [|H], |TT], [12|, [13|], it 
is also applied to pion-nucleus scattering [TJ], |T5|. In these ten years, considerable efforts 
have been devoted to extend the INC model. One of them is the inclusion of mean- 
field effects within one-body theories 



and iV-body 



16, 17, 18, 19, 20 



23, W 



dynamics |26], g7|> S Hi) e0> u3 U3- Furthermore, it becomes possible to describe the 
time-development of totally anti-symmetrized wave functions |33|, [35| . These extended 
versions of the INC model are called microscopic transport models or microscopic 
simulations, and they have been successful in studying various aspects of heavy-ion 
collisions as well as light-ion induced reactions — particle production |5B|, |S7|, [21], [B| , 
single-particle spectra p3| , collective flow |TJ| p5| , and fragment production [^7], [K], 



This paper is organized as follows: We describe our INC model in Sec. 0, and 
present a RIA model in Sec. || In Sec. f|, elementary cross sections are presented. In 
Sec. |5|, we analyze the experimental data at KEK-E176 || and discuss the reaction 
mechanism of the (K~, K + ) reactions. Finally, Section |6] contains our conclusion. 



2 Intranuclear Cascade (INC) model 

The Intranuclear cascade (INC) model we adopt here describes the propagation of 
the leading particle (the most energetic mesons at each time step). The momentum 
distribution of a nucleon in the target is assumed to be that of the Fermi gas, and the 
probability and the point of the collision between leading particle and the nucleon are 
determined by the mean free path of the propagating leading particle. This prescription 
is not fully microscopic, since the dynamical evolution of other particles is ignored. 
For example, we ignore the change of the nucleon momentum distribution and the 
secondary collision of a produced hyperon. However, this model is enough to describe 
the K + meson production, because other particles than the leading particle cannot 
produce the K + meson energetically and distortion of the target density is expected 
to be only important after the K + meson go out from the target. 

Most of cascade type models are realized by the Monte-Cairo procedure. There 
are several different prescriptions for a cascading of the colliding particles. One method 
is a so-called "closest distance approach" |12|, [l3|; if the minimum relative distance r min 



for any pair of particles becomes less than y <j(^)/tc, that pair is assumed to collide 
with each other, where u{\fs) is the total cross section for the pair at the cm. energy 
yfs. In another one, the collision frequency is determined by the mean free path. For 
example, the collision will take place with the probability a(^)vp{r)dt |10|, [H], |T5|| , 
where v and p(r) are the relative velocity of the colliding pair and the nuclear density, 
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respectively. We have checked that the both methods give the same results for the 
K + momentum distribution in the (K~ , K + ) reaction However, other details may be 
different in several computer codes, therefore we explain our Monte-Cairo procedure 
for completeness as follows: 

1. The incoming K~ meson which is the first leading particle is boosted at the 
incident momentum px- = 1.65GeV/c in the laboratory frame. 

2. To determine where this particle collides with a nucleon in the target nucleus A Z, 
we generate the Fermi momentum randomly with pp = (37r 2 pi) 1//3 , where p\ = 
pZ/A, p = Z/(4:/3iiR 2 ) and R = 1.2A 1 / 3 (fm) in the case of matter distribution for 
the nucleus. If we employ the Fermi type nuclear density, the local Tomas-Fermi 
approximation is used to obtain the Fermi momentum Pf(t). 

3. The total cross section between the leading particle and the proton (neutron) 
cr p ((j n ) is determined at this energy. If the leading particle is a resonance, the 
decay probability per unit length is calculated as Pd = r/(vy), where r is the 
total width of the resonance, v and 7 denote the velocity and Lorentz gamma 
factor, respectively. 

4. The mean free path is obtained as A = 1/(P P + P n + Pd), where P p = cr p p p , P n = 
a n p Q . We determine randomly the collision point according to the probability 
distribution dP/Ar = exp(— Ar/X)/X: This distribution is realized by choosing 
Ar as follows. 

r' = r + Ar-^- (1) 

IpI 

with p is the momentum of the leading particle, where Ar = — Alog(l — x) and 
a; is a random number between < x < 1. 

5. If |r'| is larger than the nuclear radius R, cascading is regarded as finished. 
Otherwise, we select the branch randomly according to the probabilities of P p , 
P n and Pd- 

6. If the leading particle is determined to collide with other nucleon, the collision 
branch is chosen randomly according to the partial cross sections, i.e. the elastic 
or inelastic collision between the leading particle and the nucleon. After choosing 
the branch, the kinematics of the collision is calculated, and the final momenta 
of particles after this collision are determined. 

7. If particles after the collision are nucleons, the Pauli blocking is checked. If its 
momentum in the laboratory frame is smaller than the Fermi momentum 270 
MeV/c, this collision is considered to be Pauli blocked and nothing occurs. 

8. The procedures from 2 to 7 are repeated until this event finishes. 

It is useful to mention here that our cascade procedure can be used only for the 
matter (constant) nuclear density p = po6(R — |r|). If more realistic density is used in 
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the calculations, as the Woods-Saxon type density, the collision prescription has to be 



modified. In that case, we change the collision rate according to Ref. |10|, [15] or "closest 
distance approach". We have checked that the calculated K + momentum distribution 
does not depend on the cascading prescription. 

Since the total cross sections are very small in the (K~ , K + ) reaction, note that 
a large amount of computational time is required to get sufficient statistics within 
the Monte Cairo technique. In this work, therefore, while all the elastic scatterings 
and the relevant inelastic reactions described in SecIJ are included, we ignore other 
non-interesting inelastic processes which do not result in the K + production; the prob- 
abilities for relevant inelastic reactions are increased by the factor of 

where denotes the partial cross section for relevant inelastic reactions, and a to t 
and a e i are the total and elastic cross sections at each collision energy, respectively. 
This procedure modifies the weight for each simulation, since we only treat the elastic 
scattering and relevant reactions. Thus when the relevant inelastic branch is chosen, 
the weight for the j-th simulation event after z-th step is defined as 

w\ = , (3) 

where the initial value is w J = 1. Under this assumption, the probability to detect the 
K + meson with a given impact parameter b is calculated as 

d 2 P(b) = gj w*& 1 
dpdVt N ev dVt^dp 

with dfld = 27i(cos9 min — cos9 max ), where 9 min = 1.7° and 6 max = 13.6° are the 
minimum and the maximum angles detected by the experimental measurement ||, 
respectively, and N ev represents the number of events generated in the simulation. The 
weight is defined as the final weight of the j-th event, and B J ' = 1 if the K + particle 
enters the detector range, otherwise J = 0. Typically, we generate 1,000,000 events 
in this calculations. 

In order to compare the theoretical results with the experimental data ||, we 
integrate the probability over the impact parameter b: 

i*°\'-t dh im, (5) 



\ dpdfl I J dpdVL 

where the definition of () is the forward cross section averaged over the detected solid 
angle in the laboratory frame. 



3 Relativistic impulse approximation (RIA) model 

In this section, we explain the relativistic impulse approximation (RIA) whose 
results are compared with those of INC. The RIA calculation in this paper is the 
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standard one in the quasifree region, and is essentially the same as Iijima et al. ||. 
The one-step cross section K~ + p — > K + + is calculated from the invariant 

matrix element M. and the momentum distribution of a nucleon in the target nucleus: 

da = Z eff (K~,K + ) /^^^/(p 2 ) (|-M| 2 ) S (2^)V(p 1+ p 2 -p3-P4)(6) 



d 3 Pi . ( AttP j 



where the suffixes 1, 2, 3 and 4 represent K~,p,K + and E~, respectively. As for 
the proton momentum distribution in nuclei f(j>2), we have adopted Fermi gas model 
with P F = 0.27 GeV/c. In this calculation, we use the parameterized elementary 
differential cross sections, as we will express in Sec. £|, so as to calculate the spin- 
averaged invariant matrix element squared ( | Ai | 2 ) s for elementary process of K~ +p — > 
K + + S~/S*~(1530). The elementary cross section is written as 



da elem \s) 1 1 Pf 



(\M\ 2 ) 



dttcM ^E x E 2 v X2 (2tt) 2 4^ 



s 



where yfs andp/ are the invariant mass and the final cm. momentum of the elementary 
process, respectively. Using Eqs. © (0) (f§), we get the one-step quasifree cross sections 
as follows, 

d 2 a , pi f d 3 p 2 d 3 p 4 y/s da eler " 



ry (ir - ts+\ Ps [ d°p 2 d°p i y/s da~[s) 4 

dp a Ze " {K > k) e 3 J f{P2) J; Utt^r 6 (pi+pz-ps-p*) • ( 9 ) 

In the RIA model, distortions of the initial K~ and final K + mesons are usually treated 
in the eikonal approximation. Thus the effective proton number Z e ff(K~ , K + ) which 
represents the collision frequency of the particles, is given by 

Z eff(K~,K + ) = —J rfrp(r)exp -a K - J p(x, y, z')dz' - a K + J p(x,y,z')dz' , 

(10) 

where p(r) denotes the nuclear density, and ax- = 29.0mb and &k+ = 18.4mb are K 
-nucleon and K + -nucleon total cross sections, respectively. 

For productions of the intermediate vector and scalar mesons and their decay, 
we use following expressions in RIA: 

da = Z eff (K-,M) J ^^^^ f(P2) (\M\ 2 )s (2nhr6\p 1+P2 -p 3 -p,) 

x dr 5 dr 6 ( 2 7ih)W(p ?j - Ps - P6 ) (ii) 

M = M 12 ^ u G 3 (u 3 ,p 3 ) M 3 _ 56 (12) 

where the four momentum of the meson is represented asp3 = (u 3 , p 3 ), and the suffixes 
1, 2, 3, 4, 5 and 6 represent K~,p, M, A, K + and K~, respectively. In this reaction, 
we have assumed that (1) the lifetime of the meson is long enough, and (2) the decay 
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occurs isotropically. The assumption (1) is realized by the approximation to the Green 
function G3. In the case of the scalar meson, the Green function is given by 



|G 3 (^3,P3)r 2 = |^ 2 - El + zM 3 r 3 \ 2 ~ 2 * 5(uJ 3 - E 3 ) , (13) 

3 33 



E 3 = y/Ml + pl, (14) 

where r 3 is the total decay width of the meson. Note that the assumption (1) is already 
used in Eq. (|TTD implicitly since the effective proton number Z e ff(K~, M) is calculated 
under the assumption that the meson M does not decay inside the nucleus. The matrix 
element can be calculated under the assumption (2), which is obtained as 

|AW = 8?rM ^ 56 , Vf = J M* - 4MI/2 . (15) 

Thus, the differential cross section for the K + particle through the meson decay can 
be expressed as 

d 3 a f Mm \ 2 D fdQ M ^ d 3 a(K~A — » MA) 



Br(M ^ K + K~) = r M ^ K+K -/r M , (17) 

where the solid angle Qm is calculated in the rest frame of K + , and E 3 d 3 a(K~ A — > 
MA)/dp 3 M denotes the invariant cross section for the one-step intermediate meson 
production per a proton, which is given by Eq.([5]) without the effective proton number 

Z eff- 



4 Elementary cross sections 

Now let us explain the elementary cross sections used in this paper. All the energies and 
cross sections are give in GeV and mb unit, respectively. We have categorized the events 
into three types. The first category is the direct-type reaction, which concerns the one- 
step double strangeness exchange reaction (KN —>■ KE, KB*). The second one is the 
heavy scalar and vector meson productions and their decay (^A^ — > MY, M — ► K~K + ) 
proposed by Gobbi et al. 0. The third type is the two-step strangeness exchange and 
production reaction, such as KN — > ttY followed by ttN — > KY . We will explain the 
necessary elementary cross sections in this order. 

The following processes are necessary to describe the direct-type double strangeness 
exchange reaction. 

(A) All the elastic scatterings and charge exchange reactions of KN and A^A^. 

(B) Double strangeness exchange processes: 

KN -> KE, ATH*(1530). 
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(C) Initial anti-kaon interaction: 

KN -> K*(892)N. 

(D) Decay of resonance kaons (anti-kaons): 

K*(892) -> fGr, A"* (892) -> Ati, 

(E) Final kaon interactions: 

KN — > (A*(892)iV, AA(1232)), 
A*(892)iV -»■ (if A, A*(892)A(1232)). 

Although the actual strangeness exchange occurs only in the processes of (B), initial and 
final state interactions have to be included when we are interested in the absolute and 
whole spectrum of K + including the lower momentum region, as will be demonstrated 
in Sec. [|. 

The above cross sections are fitted to the experimental data [R^, [4(J with appro- 
priate functions or tabulated. The total, elastic and charge exchange cross sections of 
(A) are well-known experimentally [[39], f§T | . We employ here the smooth interpolations 
of these data. Angular dependencies of the elastic collisions are simply assumed to 



have the form 12, 16 



da 
~dt 



oc exp(b(y/s)t) 
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where t denotes the square of the transferred four momentum, y/s is the cm. energy, 
and b(yfs) is taken to be b(y/s) = 0.125a (y/s) to t- 

The double strangeness exchange reactions (B) are also well-known experimen- 
tally although the errorbars are still large. Therefore, we use following parameterization 
for the elementary H and S*(1530) production cross sections: 

0.0413463(^i - 1.81496) 168792 



<T(K~p -> K + E~) 



(y/s- 1.93536) 2 + 0.01174 
0.00613(v^- 1.8 1 496) - 48020 

(y/s- 2.02435) 2 + 0.01739 
0.001673(v^ - 2.11496) ' 338065 



[mb] 



v^< 2.22[GeV/c] 



m 



b] 



a(K~p -> K U E , 
a(K~p -> A+H*-(1530)) 
a(K~p -> A H*°(1530)) 



2.3174) 2 + 0.016429 

0.00364308^- 1.8 1 25 71) - 620511 

(y/s - 2.05469) 2 + 0.0184015 
0.0017692(v^ - 2.0264) - 45434 



y/s> 2.22[GeV/c] 
(19) 



[mb] , (20) 



(^- 2.09274)2 + 0.01394 
0.003854(v/i - 2.0264) 3475 



[mb] , 
mb] . 



(21) 



(22) 



(y/s- 2.18855) 2 + 0.028487 

These fitted cross sections are shown in Fig. [L], together with the experimental data. 
Although we can take another parameterizations within the experimental error bars, 
we have checked that our simulation results are almost the same. 
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Fig- 1 

Around the incident momentum of px- = 1-65 GeV/c for the K~ beam, the experi- 
mental data shows a backward peak angular distribution. We fit angular distributions 
of S production cross sections to the experimental data at p K - = 1.7 GeV/c, neglecting 
these energy dependence: 

At f 0.34308 exp(0.6229 cos 6) + 0.06681 exp(-4.4396 cos 6) for K~p -> #+5" , 
dfi ~ \ 0.44 exp(3.1684 cos 6) + 0.5428 exp(-1.6747 cos 6) for K~p -> . 

(23) 

As shown in Fig. |2|, this parameterization is a good approximation within a Fermi 
momentum spreading around px- = 1.65 GeV/c. For the S*(1350) production, the 
angular distributions are assumed to be isotropic. 

Fig-j 

In order to describe the reactions in the second category, following elementary 
cross sections and the decay widths are required. 



(F) Scalar and vector meson production: 

KN^(<f>,a J )A. 

(G) Decay of scalar and vector mesons: 

^,a ,fo)^K~K+. 

Compared to the processes of (A)-(E), the experimental information on the scalar 
and vector meson (</>//o/ao) production is very scarce. Here, <fi/ fo/ao production cross 
sections and angular distributions are taken from Ref. 0: 

a(K~p — >■ + A) = 0.31531p/exp(-1.45p/) [mb], (24) 
a{K~p -> f + A) = 1.51a(K-p-+ <p + A), (25) 
a(R-p^ a + A) = 1.38a(K~p —></> + A), (26) 

exp(1.8cos0) , (27) 



da 

In 



where pf is the relative momentum for the final state. The angular distribution of the 
(f> decay into KK is assumed to be isotropic in its rest frame. But we have checked that 
an anisotropic decay (cos 2 9 distribution) has no significant influence in our results. 

For two-step processes, in the case of the 12 C target, Iijima et al. found that 
the calculated cross section is smaller than the observed one in an order of magnitude. 
Moreover, Gobbi et al. || also suggested that two-step processes can not reproduce 
the experimental cross section at low momentum region in the K + spectrum. Note, 
however, that the intermediate meson is limited to the pion, and they considered only 
for the 12 C target case. We expect that inclusion of heavy mesons as rj mesons would 
strongly affect the final results especially in the case of heavier nucleus. Therefore, we 
take account of two-step processes on all targets by including following processes: 
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(H) Single strangeness exchange reactions: 



KN^(ir,p,ri,u,rf)(Y,Y*) . 

(I) Single strangeness creation reactions: 

(ir,p, V ,u, V ')N^(K,K*)(YY*), <j>N . 

(J) Decay of resonances: 
p — > 2tt, uj — > 3tt. 

Here, Y = {A,£}, Y* = {A(1405), A(1520), £(1385)} and K* = {K*(892)}. We treat 
explicitly the isospin of all these particles in our calculations. 

We have fitted the experimental data as far as they are available. The fitted 



function is partly taken from Refs. fl50| , [52jj (see Appendix). Other unknown cross 



sections in the second step (I), are estimated using the Breit-Wigner formula in the 
same way as Ref. || ||. Namely, processes of (I) are assumed to occur through A^*, A* 
resonances and interferences of various resonances are ignored. Since the energy is 
located just at the baryonic resonance region in the {K~, K + ) reaction with the incident 
momentum of p%- = 1.65 GeV/c, this assumption is expected to be valid. The explicit 
form is expressed as 

a(MB _> M'B') = JL V (2^ + 1) r R (MB)r R (M'B') 
1 1 kl m Z ^{2S M + l){2S B + l){^S-m R Y + r R {totY/ A 1 1 



with the momentum dependent width 

r R (MB) = ( A)' 2 '"' f r&MB) , 

\p R J l + 0.2(p/p R ) M 
where r^(MB) denotes the partial width of the R — > MB decay, and 



(29) 



P = 7T-r= , (30) 





(m B + m M ))Ws- 


(m B - 


m M )) 






(m B + m M ))(m R - 


- {m B - 


- m M )) 



PR = ^ , (31) 

are the relative momenta between the meson and the baryon in the final state at a given 
cm. energy and the resonance mass, respectively. The summation R in Eq. (p8|) runs 
over resonances, as N(1440) ~ JV(2190), A(1600) ~ A(1950). This formula is called 
as the generalized Breit-Wigner formula, which is obtained by neglecting t channel 
contributions and the interference of the resonances. Actual values for these parameters 
are taken from the Particle Data Group ||41|| , as listed in Tables [l] and 0. 

Table. | 



10 



Table. | 



The formula Eq. enables us to estimate experimentally inaccessible cross 
sections such as pN — > AK. The N* resonances decaying into the uN and r/'N branches 
are included to fit the experimental data of ir~p — > ujn(i]'n) |3H[| . It has been shown 
that such meson resonances play an important role to product strangeness particles in 
heavy-ion collisions at AGS and SPS energies |5j. In fact, strangeness production cross 
sections via the meson resonances are found to be much larger than ttN cross sections. 

Although t channel component of KN interactions are considered to be large 
in (H), we also apply the resonance formula of Eq.(|28|) taking into account hyperon 
resonances A(1520) ~ A(2350) and £(1660) ~ £(2250) in order to obtain the cross 
sections. However, there is a lack of experimental information about resonance pa- 
rameters compared with N* and A*. Therefore, we determined the experimentally 
unknown resonance parameters to fit the existing data such as K~ p — > 77 A. The 
parameters are shown in Tables |3| and f|. This procedure would be sufficient for this 
study because we need only KN interaction cross sections as the incoming states. If 
we cannot obtain the cross sections and angular dependencies from the Breit-Wigner 
formula or from the experimental data, for example, cr(p~p — > £°Af°(892)), we simply 
assume that these values are the same as ttN cross sections at the same cm. energies. 
For the total cross section of 4>N, we have used the value of 8.3mb, which was suggested 
in a photoproduction experiment P3"| , while 16.0mb is predicted in the additive-quark 
model Ref. @. 



Table. g| 
Table. | 

Note that the process 4>N — > KA is also able to produce the K + meson. The 
cross section calculated from the OBE model of Ref. amounts to about 5mb, as 
seen in Fig. [3|. Our INC calculation includes this process as a two-step process. 

Tig~|" 



5 Results and Discussions 

5.1 S and H*(1535) Productions 

Let us now compare our results with the experimental data. At first, we focus on the 
S and S*(1535) production processes. In Fig. |], we display the INC spectra using both 
the matter density (solid lines) and the Fermi-type nuclear density (dashed lines) for 
(K~, K + ) reactions on 12 C, 27 A1, 63 Cu, 107 Ag and 208 Pb targets. The RIA results with 
matter density are also presented (dotted lines). 
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Fig. 1 

We can see that the results with the Fermi-type density well reproduce the exper- 
imental data at the high momentum region, while the results with the matter density 
underestimate the experimental data. Our calculations with Fermi-type nuclear den- 
sity for 5 and S* production processes are consistent with DWIA of Ref. ||, which 
was calculated by using the same nuclear density jKJ. It is also seen that the INC 
and the RIA results with matter density give almost the same results around the high 
momentum region of K + , including the shape and absolute value of the spectrum. 
The difference can be found in the lower momentum region of K + . This is due to the 
different treatment of initial and final state interactions. The eikonal approximation of 
Eq. ([[(]) means that once the particle is scattered, the flux is assumed to be lost in the 
RIA model. On the other side, the INC model treats the initial and final state interac- 
tions explicitly, which ensure us that both loss and gain of the flux are automatically 
included. For example, some parts of the scattered K + particles will be detected while 
the momentum becomes lower than that just after the production. 



5.2 Contributions of the scalar and vector meson decays 

Here we consider the effect through the scalar (a /f ) and the vector meson((/>) decays. 
The importance of this process were first suggested by Gobbi et al. ||. In Fig. [5], we 
present the calculated results of the K + spectra (solid lines), taking into account the 
decay from (p, a and f mesons; the (ft/ fo/a production and their decay into K~K + 
significantly contribute to the K + spectrum at low momentum region. However, these 
processes are not enough to reproduce the large bump in the experimental data. 



Fig- | 

In order to test these results, we also calculate the cross section of the above 
processes employing RIA. The RIA results are shown in Fig. |5| (dotted lines). We find 
that INC and RIA calculations give us similar results for all targets. However, our 
results are about 50% smaller than the cross section in Ref. |J. There are several 
probable reasons for this contradiction. The one reason lies in the treatment of the 
Fermi motion in nuclei. Since the (ft/fo/^o mesons can be produced by a subthreshold 
particle production in the (K~,K + ) reaction at p K - = 1.65 GeV/c, the effect of the 
Fermi motion is expected to be very important in this situations [|7[ |2T|. Thus we have 
checked the sensitivity of the treatment for the Fermi motion. As seen in Fig. |6|, we 
compare the results of following calculations for the 27 Al target with the matter density 
(matter), the local Fermi momentum with the Fermi type nuclear density (F+F) and 
the matter Fermi momentum with the Fermi type nuclear density (F+M). It can be 
seen that the cross section in case (F+M) is enhanced by the factor 1.3. This treatment 
corresponds to the calculation in Ref. ||. 

"FiglT 
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In addition, we mention that Gobbi et al. || used an "effective" incident mo- 
mentum which includes the effects of nucleon Fermi motion, and regarded it as the 
incident momentum. If this effective incident momentum is used in the original lab- 
oratory frame, it gives a wrong kinematics. In fact, the high momentum threshold of 
the K + particle is artificially modified to be about 1.1 GeV/c in Ref. 0, while the 
actual threshold is about 0.8 GeV/c through 4>/a / f decays. 

5.3 Two-Step Processes 

In Figs. [7| and ||, we show the calculated K + spectra through two-step processes 
(dotted lines) on 12 C, 27 A1, 63 Cu, 107 Ag and 208 Pb targets: Long dashed lines denote 
the contributions of the 5 and S* production, dashed lines represent those from (p/a / f 
decay, and solid lines are the sum of all the processes considered in this paper. It is 
clearly observed that two-step processes play an significant role in K + yields, especially 
for heavier targets. This result is in contrast to Ref. [J. In addition, the calculated 
total spectra well reproduce the experimental data. 



Fig. 



Before concluding that two-step processes are important, we should estimate the- 
oretical ambiguities in the present calculation. We consider that the largest ambiguities 
come from the angular distributions of elementary processes. Here we have used the ex- 
perimental angular distributions if possible, for example, the S, S*, productions and 
some part of two-step processes. The experimentally unknown angular distributions 
of two-step processes such as pN are simply assumed to be the same dependence to 
that of tiN at same c.m.energy. Therefore, it is necessary to estimate the sensitivity of 
the angular distributions. In Fig. |fj, we represent the results of the two-step contribu- 
tions assuming the isotropic angular distribution for all of the two-step reactions. This 
choice decreases forward angle cross sections, but the difference of integrated yields is 
at most 30%, and their contributions are still significant. Due to the large uncertainty 
of the elementary angular distribution like pN — > KY, more quantitative discussions 
are impossible at this stage. However, we can conclude that a significant part of K + 
particles are produced though various two-step processes, since the use of isotropic 
angular distribution might give the minimum yields of K + . 



Figjj 

Contrast to the previous works || the reasons why two-step processes have 
large cross sections can be understood as follows. The first reason is that we take 
account of the meson resonances p, i], uj as a intermediate meson. In addition to the 
increase of the number of processes resulting in K + production, due to the large masses 
of meson resonances, the strangeness productions through these mesons have a larger 
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cross section than that through pions. Note that the irN reaction for strangeness 
productions is all endoergic one (Q < 0), whereas p-, 77- and cu-nucleon reactions for 
strangeness productions are exoergic one (Q > 0). Thus inclusion of these mesons is 
of importance. Figure [10] displays the contributions of intermediate mesons. It can be 
seen that tt, 77 and u mesons mainly contribute to the two-step cross sections. It is of 
importance to see that the cross sections of meson resonances are comparable to that 
of pions. The reason of the small p-contribution is that the p meson are likely to decay 
before its collision due to its large width. We find also that pions from the p meson 
decay can contribute to yields of K + . We confirm that if we include only pions as 
the intermediate meson and neglect hyperon resonances, this calculation gives similar 
result to the estimation of Iijima et al. in the case of 12 C target: As expected from 
the small mass number, the contribution from the two-step processes is rather smaller 
than that for the heavier targets. 



Fig. 10 



The second reason is related to the meson momentum region. The mesons pro- 
duced by K~ at p K - = 1.65GeV/c are just in the baryon resonance region where the 
cross sections have the largest values. Alternatively, at the higher momentum of K~, 
the contribution of the <fi meson and its decay becomes larger, and K + from the (f> 
meson decay would dominate K + yields. 



We can also consider the 4>N — > K + A process. Using the OBE model f4"l 



we 



can obtain the cross section about 5mb at most. Although this process is taken into 
account in the INC calculation, it is found that the contribution is small. 



6 Summary and Conclusions 

In this paper, we have analyzed (K~, K + ) reactions within a INC model that is newly 
developed especially to include various inelastic channels in order to investigate the 
whole momentum region of the K + momentum. The detailed comparison for the 
K + spectrum has been done between the INC and the RIA model calculations. We 
have shown that the results with these two models reasonably agree with each other 
in one-step processes, and the calculated total spectra with INC well reproduce the 
experimental data on various targets. 

Our results show that the momentum distribution of K + can be explained mainly 
by following mechanisms: 

(1) Direct-type reactions: 

K-p -> K+E- , 
K-p-^ A^+H*-(1530) . 

(2) Decay of scalar/vector mesons: 
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fo J 




{ fo 



K-p - a A, a ^ 



(3) Two-step processes: 



K-N followed by MN - | ^ | j £ 

where M = n, p, i], oj. 



For the process (1) and the process (2), the consistent results are obtained between the 
INC and the RIA models. Among two-step processes (3), it is important to include 
it, p, t], uj mesons, and hyperon and kaon resonances to reproduce the K + spectrum. 
There are mainly two reasons of the enhancement of the two-step contributions, con- 
trary to the estimation in the previous works |], The first one is the variety of 
intermediate mesons and final hyperons. It is clear that a number of the path resulting 
in K + production is much larger, because various mesons and hyperon resonances are 
included. In addition, meson resonances have larger masses than pions, and stored 
energies in masses are released in second step reactions. Thus the K + production cross 
section with these meson resonances becomes larger than that of pions. The second 
reason is related to the meson momentum region. The mesons produced in this reac- 
tion is just in the baryon resonance region where the cross sections have the largest 
values. 

In conclusion, the INC model can explain consistently the K + momentum spec- 
trum from the (K~,K + ) reaction at p^- = 1.65GeV/c on various targets. Both 
the two-step strangeness exchange and production processes with various intermedi- 
ate mesons and 0, ao and fo productions and their decay into K + K~ of which was 
first quoted in Ref. |J, are necessary to reproduce the experimental data in the low 
momentum region. 

There are still some theoretical ambiguities and problems in the present study. 
One of them is, as already discussed in Sec. [5], the angular distribution of elementary 
processes. Another problem is that we have ignored the propagation of heavy baryon 
resonances which are assumed to be formed in Eq. (^); however, the lifetimes of these 
heavy resonances are so small that they are expected to decay before interacting with 
other nucleons at around the normal nuclear density. The third one is the branching 
ratios of heavy baryon resonances, especially decaying into uN. Since the measurement 
of decaying into uoN is relatively difficult, the experimental information of the branching 
ratio is scarce. Although we have fitted the u production cross sections and determined 
the plausible branching ratios, it may be problematic to assume this reaction occurs 
only in s channel. These facts may be related to the overestimation of the K + spectra 
at very low momentum region. 

There are some interesting problems suggested in this work. First, we note that 
the (K~ , K + ) reaction can produce hyperon resonances Y* in the nuclear medium. This 
opens up the possibility of studying the properties of the hyperon resonances in nuclear 
matter. Another interesting problem is to predict yields of double-A hyperfragment 
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production and to make their production mechanism clear. This is possible if we 



extend our cascade code including mean field effects p(| £7], £3, |T|, |3(J, |31], |32|. In 
order to search for double-A hypernucleus experimentally, the (K~,K + ) experiment 
is planned at AGS(E906) [f48| , [49[| in the near future. The production mechanism of 
hyperfragments may be different from the K + productions through direct S, direct S*, 
decay and two-step processes. The study in this line is now in progress. 
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A Parameterizations of Cross Sections and Angu- 
lar Distributions 

A.l Cross Sections 

We list up the cross sections parametrzied with the experimental data. Some parame- 



terizations are taken from Ref. [50, 51, 52 



(K~ : A ') °-° 2(v/i ~ ^ )0A 0.000186093^ - ^gg 

a{ V) (^- 2.20726)2 + 0.0589288 (y/s - 3.07199) 2 + 0.00652879 1 ' 

_ 0.007665( v ^-1.613) a1341 

a 7r p — » AK ) = , — — 33 

V v ' (y/s- 1.72) 2 + 0.007826 V ; 

^ r,+s 0.009803(7?- 1.688) - 6021 0.006521^ - 1.688) 1 - 472 * A . 

a(7i~p -> T,~K + ) = vv ; _ — — H vv , _ — -f34) 

v y ' 0.006583 + (y/s- 1.742)2 0.006248 + (y/s - 1. 940) 2 " ; 

,(0> - 2°A-») = 0-05014(>A -l.«88)'-»" 



0.006455 + (y/s - 1.730) 

Xr ,„ N 0.02655(v^-2.018) a5378 , , 

,(,-,, ~A(1405)tt«) = UT -^ fwT J m (36) 

.(,- p A(1520)Jf) = 0W«763(Vi - 1.833)""- 

v y v ; ; (^-2.21759)2 + 0.0249613 v ; 

a(7r-p -> AAT°(892)) = 0.240985gexp(-1.09397g), (38) 



'(a - (M A + 0.636))(a - (M A - 0.636)) 
where q = -= (39) 

2y/ S 
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x , N 0.00163105(Vs - 1.752) , x 

a(,-p -» ^-(1885)^) = (v - s _ 2 1833 y + Q on49 ' 82 (40) 

„ n/ xr , nx 0.0113016(^- 1.752) 
,(0, S°(1385)^) = (v - _ 2 124 ^ 2 + Q (41) 

a(n-p -> £-#+(892)) = 0.123117gexp(-1.50466g), (42) 

^ - (Af E - + 0.636))(s - (Me- - 0.636)) 

where q = -= (43) 

2y s 

a(n-p -> S°ir°(892)) = 0.188256gexp(-1.32449g), (44) 

a/(s - (M E o + 0.636)) (s - (M s o - 0.636)) 

where g = ——= (45) 

2y s 

. _ 0.000986473(v^- 1.948) a045244 
a 7r » — > n0 = r~F= ^ (46) 

V y V! (y/s- 2.09) 2 + 0.0290906 V ; 

+ 0.03591(v^- 1.688) - 9541 0.1594(^5 - 1.688) - 01056 

a ^ P ~* >~ 0.01548 +(v^- 1.890) 2 + 0.9412 + (^ - 3.000) 2 ( ' 

*(«+p -+ £+(1385)*+) = 0-0191^(v^- 1.833)Q-~ 

V y y ' ' (y/s- 1.98968)2 + 0.0253142 v ; 

V y V ;; (v/s- 2.37123) 2 + 0.177898 V ; 
. 0.003978(v^- 1.688) - 5848 0.04709^ - 1.688) 2 - 1650 

<T(7T p — > 2j A J = - nn „„- n p 7= ^77-^ + n nnr0 r , ? 7= 1 nntr\9 ( 5U J 



0.006670+ (v/s- 1.740)2 0.006358 + (y/s - 1.905) 2 
0.05014(^- 1.688) 1 - 2878 
0.006455 + (y/s- 1.730)2 



.(A = ° n °i°! 4( ^- 68 ! i) 'Z (so 



A. 2 7riV Angular Distributions 

The angular distributions are fitted as the following function, 

^ ~ a exp(7i cos 0) + b exp(— 72 cos 0) + c(l — cos 2 9) (52) 
at the cm. energy s in GeV. 



1. 7r p 



b = 



■AK° 
0, 



-5616.87 + 6491.43s - 1859.92s 2 , 
0.0962773 

(s- 1.89039)2 + 0.00650021' 
0.0 

-155.986 + 240.621s 
-4.72881 + s 2 - 20857 , 
-65.8086 + 64.8666s 
0.0, 



s < 1.878, 
s > 1.878, 



(53) 
(54) 



123.717s 2 + 21.2028s 3 , 
15.8761s 2 , 



s < 1.878, 
1.878 < s < 2.024, 
2.024 < s < 2.097, 
2.097 < s < 2.35, 
s > 2.35, 
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j -9.54489 + 6.20974s, s < 1.878, , , 

71 7 \ 28.4659 - 33.8566s + 10.3904s 2 , s > 1.878, ^ ' 

7 2 = 83.5865 - 55.3861s + 8.68325s 2 . (56) 



2. n-p^Z~K+ 

c = 0.0 (57) 

f 0.0049, s < 1.978, , . 

a ~ { 0.0, s > 1.978, (58) 

f 0.3632, s < 1.978, 

b = { 1.0, s > 1.978, (59) 

71 = 4.32675, (60) 

1.38911, s< 1.978, 

72 = <j 378.699 - 361.502s + 86.4954s 2 , 1.978 < s < 2.097, (61) 
-8.44056 + 4.53549s s > 2.097. 

(62) 



oz^o 



c = 0.0, (63) 

a = 0.00537342, (64) 

6 = 0.391835, (65) 

7! = 5.17482, (66) 

72 = 0.582968. (67) 



4. 7T+p -> 



0.00806605, s < 1.9, 

0.000323596, s > 1.9, 

10.0601, s < 1.9, 

0.00018645, s > 1.9, 



(68) 
(69) 



f 0.0, s < 1.9 , . 

; : \ 38.3724 s > 1.9, [ > 

} 7.49409, s < 1.9, 

ri \ 13.3355, s > 1.9, 1 ; 

f 0.889499, s < 1.9, , . 

72 ~ \ 13.3355. s > 1.9. ^ ^ 

(73) 
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Table 1: Branching ratios for A* decay 





J p 


Mass 


Width 


N*7C 


Ntt 


Avr 


Np 


UK 


N* type 


-P33 


3 + 
2 


1.60 


0.35 


0.15 


0.175 


0.55 


0.125 


0.00 


N(UA0) 


5*31 


1- 
2 


1.62 


0.15 


0.05 


0.25 


0.50 


0.20 


0.00 


iV(1440) 


£33 


3- 
2 


1.70 


0.30 


0.00 


0.20 


0.40 


0.3984 


0.0016 




5*35 


1- 

2 


1.90 


0.20 


0.40 


0.10 


0.05 


0.45 


0.00 


iV(1440) 


-^35 


5 + 
2 


1.905 


0.35 


0.00 


0.10 


0.14998 


0.75 


0.0002 
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0.25 


0.60 
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0.25 


0.15 


0.00 
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0.03 
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0.01 
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Table 2: Branching ratios for N* decay 
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£>15 


5~ 
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0.34 
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Table 3: Branching ratios for A* decay 
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Table 4: Branching ratios for E* decay 
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Figure 1: Parameterizations of (a)S and (b)5* production cross sections with experi- 
mental data [02 . 
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Figure 2: Angular distributions of (a) K~p — > K + E~ and (b) K~p — > K°E°. The 
solid curve denotes for the angular distribution parameterized at the laboratory K~ 
momentum of pk- = 1.7GeV/c. 
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<|> N -> KA cross section 
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Figure 4: Comparison of the calculated results for the one-step processes (S~ and S*~ 
production) with the experimental data. Solid and dashed lines show the INC results 
with the matter density and with the Fermi-type density + local Fermi momentum, 
respectively. Dotted lines show the RIA results with the matter density. 
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Figure 5: Comparison of the calculated results for the meson-decay(0 + ao + /o) contri- 
butions to K + momentum spectra with the experimental data. Solid and dashed lines 
show the INC results with the matter density and the Fermi-type density, respectively. 
Dotted lines show the RIA results with the matter density. 
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Figure 6: (K~,K + ) double differential cross sections on 27 A1 target at p K + = 1.65 
GeV/c. The INC calculations only contain 0, ciq and /o productions. These spectra 
are calculated by using matter density (solid line), the Fermi density with the local 
Fermi momentum (F+F) (long dashed line) and the Fermi density with the matter 
momentum (F+M) (dashed line). 
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Figure 7: Calculated momentum spectra of K + for 12 C, 27 A1, 63 Cu, 107 Ag, 208 Pb and 
proton targets at px- = 1-65 GeV/c using the INC model. The matter density is 
used. The squares represent the data of Iijima et al. [[J. The contributions of H and 
5(1535) productions are represented by long dashed lines. The dashed and dotted lines 
correspond to the contributions of 0, a and fo productions and of two-step processes, 
respectively Solid lines denote the results of the total spectrum. 
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Figure 8: Calculated momentum spectra of K + by using the INC model for various 
targets at px- = 1-65 GeV/c with the Fermi type density for targets. The notation is 
the same as fig [7|. 



31 




0.0 0.5 1.0 ^ 1.5 

K + momentum (GeV/c ) 



Figure 9: (K~, K + ) double differential cross sections for only the two-step processes in 
the INC model, together with the experimental data of Ref. ||. Solid lines correspond 
to the results with isotropic angular distribution for the MN — > K(K*)Y(Y*), M = 
7r,p,r],uj and rj . Dashed lines represent the results in which resonance- N angular 
distributions are equated to those of irN scattering. 
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Figure 10: Decomposition of the two-step processes calculated with the INC model. 
The contribution from ir, p, 77 and u intermediate mesons are represented by solid, long 
dashed, dashed and dashed-dotted lines, respectively. 
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